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Abstract 

Estimates for leading and non- leading 'twist' distribution functions are obtained within the 
framework of a diquark spectator model using a non-local operator representation. 
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More details about the method and the results reported here can be found in the long write-up [T|]. 

On ■ 

3£ 1 introduction 

o . 

r*- • In a field-theoretic description of hard scattering processes the information on the hadronic 



structure is contained in matrix elements of non-local operators (parton correlation functions). 
For instance, the distribution functions of quarks in hadrons are obtained from these hadronic 
V^ ' matrix elements by tracing them with certain Dirac matrices and integration over components 
q_) | of the quark momentum. 

^ 1 Distribution functions are universal in the sense, that they occur in the same form in the descrip- 
^ ■ tion of all hard processes involving the same kind of soft physics. Being, however, of genuine 
^ . non-perturbative nature, in general, they can not be calculated by perturbative means (with 
^ ! the exception of heavy quark distribution functions). Ultimately, one could wait for lattice 
gauge theory to provide the answer, but regarding the enormous difficulties those techniques 
encounter, in the meanwhile, model estimates for the functions may be very useful in predicting 
cross sections or asymmetries in future experiments. 



2 correlations and distribution functions 

quark-quark correlation function 

The simplest correlation function containing information on the distribution of a quark (with 
momentum p) inside a hadron (characterized by momentum P and spin vector S) can be written 
for each flavor as [@, |[ §] 

*(p;P,S) 
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^(p } P,S) = ^J (P,S\^(0)\X)(X\Ux)\P,S). (1) 

inclusive lepton-hadron 

We encounter the quark-quark correlation function in the field theoretical description of, e.g. 
deep inelastic lepton hadron scattering, whose amplitude squared is depicted by the famous 
hand-bag diagram (the usual invariants indicated besides it) 

^ q i * 

2 2 



<r = -q 



Pit 1 1 p r>2 



Q 

X B 



2P-q 



In fact, in the cross section the correlation function occurs traced with a Dirac matrix, 7 + , and 
integrated over three components of the quark momentum 



h{x) = ^Jd p -d 2 p T Ti(^- 



(2) 

p+=xP+ 

which defines a distribution function, fi(x), depending only on the longitudinal momentum 
fraction x = p + /P + . The parton model predictions for the structure functions Fx and F 2 are 
given (to lowest order) in terms of the distribution functions as 

2 Fl (x B ) = ^^ = j:elmx B ). (3) 

X B a 



We introduce the notation 

(4) 



$[ r l( x ) = - f dp-d 2 p T Tr ($r) 

2 J P +=xP+ 

for the quark-quark correlation function traced with a certain Dirac matrix and integrated over 
dp~d 2 p T . 

leading twist distribution functions 

Depending on the Dirac matrix involved, different spin properties of the hadronic structure are 
probed. To leading order (in an 1/Q expansion) the following projections occur in the description 
of hard processes (A is the helicity of the hadron, St the transverse part of the spin vector) 



® h+] (x) = fx(x) (5) 

$[7+75]^) = Xgi (x) (6) 

$K* + 75]( X ) = S^hx(x) (7) 



which have an intuitive probabilistic interpretation.^ 
2 The notation of the functions follows Refs. || |). 
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interpretation 



/i(x) gives the probability of finding a quark with light-cone momentum fraction x in the 
"+" -direction (and any transverse momentum). 





gi(x) is a chirality distribution: in a hadron that is in a positive helicity eigenstate, it 
measures the probability of finding a right-handed quark with light-cone momentum frac- 
tion x minus the the probability of finding a left-handed quark with the same light-cone 
momentum fraction. 



§11 = 



h\{x) is a transverse spin distribution: in a transversely polarized hadron, it measures 
the probability of finding quarks with light-cone momentum fraction x polarized along 
the direction of the polarization of the hadron minus the probability of finding quarks 
oppositely polarized. 



h, = 



twist 3 functions 

The subleading ('higher twist') functions have no intuitive probabilistic interpretation. Never- 
theless, they are well defined via the quark-quark correlation function traced with appropriate 
Dirac matrices. The pre-factor Mj P + behaving like 1/Q in a hard process signals the sub-leading 
(i.e. 'twist' 3) nature of the corresponding distribution functions 




<S> [1] (x) 


M 
P+ 


e(x) 


(8) 


$[f a T*]( a; ) 


M 
P+ 


S% g T (x) 


(9) 




M 
P+ 


A h L (x) 


(10) 



transverse momentum 

In some processes additional information is gained by considering observables which depend on 
the transverse momentum of an external hadron (like differential cross sections). Typically, these 
are processses with at least three external momenta, since those cannot be all collinear, in general. 
The transverse momentum dependence of the observables is related to the transverse degrees 
of freedom of partons in the hadrons. Let us examplify the relation for one hadron-inclusive 
lepton-hadron scattering (see e.g. ||). The amplitude (squared) is depicted by 
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[q- jq + ,q T ^0] 




In a frame where the target momentum, P, and the momentum of the observed hadron in the 
outgoing channel, Ph, are collinear, the photon momentum will have a transverse component (as 
indicated besides the diagramm).[] 

The external transverse momentum, q T is related via momentum conservation, i.e., a S(q T + 
p T — kx) function at the vertex, to the transverse momentum components on the quark lines. 
Thus, observables differential in the transverse momentum will involve functions like 



fi(x,p T ) = - J dp Tr($7~ 



p+=xP+,p T 



(11) 



depending on the longitudinal momentum fraction and the transverse momentum components 
which are not integrated out. E.g., the differential cross section reveals a behavior 



da 



dxB dydz d 2 q 



r 



J d 2 p T d 2 k T 5(p T + k T - q T ) fi(x B ,p T ) D x (z,k T ) 



(12) 



where Di(z,kr) is a transverse momentum dependent fragmentation function which describes 
the hadronization of the quark. In straightforward analogy we define the functions 



$[r]( 



x,p T ) 



dp-Tr ($r) 



p+=xP+,p T 

which differs from definition (W) only by the left-out integration over transverse momentum. 



(13) 



leading twist distribution functions (p T dependent) 

Not surprisingly, there are more functions involved due to the additional degrees of freedom. To 
leading 'twist' there are six functions defined by 



^ + \x,p T ) = h(x,p 2 T ) 

&^(X, PT ) = Xg 1L ( Xj p^+g 1T ( Xj p 2 T )P^: 



^ a+] (x,p T ) = S^h 1T (x,p 2 T ) + ^ 



X hi L (x,p 2 T ) + hi T (x,p 2 T )^j^ 



(14) 
(15) 
(16) 



3 The reasoning in other frames is similar, one of the external momenta unavoidably will have a transverse 
component. 



Of course, all those Dirac projections reduce after integration over transverse momenta to the 
forms shown before, like for instance, / d 2 p T $' 7 ^(x,p T ) = $b + ](a;), etc. 



twist 3 functions (p T dependent) 

Similarly, there is a larger number of 'twist 3' functions depending on x and p T which are 
obtained by the projections: 



$> [1] {x,p T ) = 


M / 2^ 


(17) 


^ i] (x,p T ) = 




(18) 


<S> l ^ 5] {x,p T ) = 


m s T 2 






+ ^(^,j4) + % 5t ^p§.)) 


(19) 


& ia '^ 5] (x,p T ) = 


^7 + ^/4(^) 


(20) 


$ [i<T+ ~ 75] (x,p T ) = 




(21) 



3 spectator model 

The purpose of our investigation is to obtain estimates for the non-leading (i.e. 'twist 3') distri- 
bution functions which are experimentally poorly (or not at all) known at present. 
To this end we employ a rather simple spectator model with only a few parameters. After 
fixing the parameters by phenomenological constraints we check that the gross features of the 
experimentally well-known leading 'twist' distribution functions fi(x) and gi(x) are satisfactorily 
reproduced. 

3.1 ingredients of the model 

The ingredients of the model are indicated below f\. 

• The basic idea of the spectator model is to assign a definite mass to the intermediate states 
occurring in the quark-quark correlation functions 

pft . Ijp 



*(p;P,S) 



{P-pf = M R 



2 



Hot more detailed information about the model see also previous publications Jj], 0, ||, ^ 
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The quantum numbers of the intermediate state are determined by the action of the quark 
field operator on the hadronic state \P,S), i.e. they are the quantum numbers of a diquark 
system 





spin 


isospin 


scalar diquark 








axialvector diquark 


1 


1 



The matrix element appearing in the RHS of ([I]) is given by 

(X.\ih(0)\P,S) = (-r^-) Tfc U^S) (22) 
in the case of a scalar diquark, or by 

<X>,(0)|P, S) = C (-^—) T% Ut(P, S) (23) 



m 



ik 



for a vector diquark. The matrix elements consist of a nucleon-quark-diquark vertex T(N), 
the Dirac spinor for the nucleon Ui(P, S), a quark propagator for the untruncated quark 
line and a polarization vector e*^ in the case of an axial vector diquark. 

For the nucleon-quark-diquark vertex we assume the following Dirac structures []: 

i 

quark 




diquark 
(singlet) 



diquark ^ * \/3 

(triplet) 



T S (N) = 1 g s (p 2 ) (24) 



T^(A0 = ^75 (7 M + ^) (25) 



The functions g s / a (p 2 ) are form factors that take into account the composite structure of 
the nucleon and the diquark spectator ||. We use the same form factors for scalar and 
axial vector diquark: 

2 2 

9s /a (p 2 )=N ^ 2 _ Ar . (26) 

The flavor coupling of the proton wave function from a scalar diquark (S ) and an (ax- 
ialvector diquark with isospin component I 3 = or 1 (A or A x , respectively) 

\p) = -L |„ So ) + ^=\uA )-^= \d A x ), (27) 



leads to the flavor relations 



fx = \f' + \ft ( 28 ) 
ft = fi (29) 



3 a special case of the most general form given in 
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and similarly for the other functions. The coupling of the spin has already been included 
in the vertices. 

Putting all ingredients together analytic expression for the quark-quark correlation functions are 
obtained 



x(^ + m)(/+M) (l + a Rl5 $) UMm) 
from which the distribution functions can be easily projected. 

3.2 fixing the parameters 

The parameters of the model are fixed as follows: the power in the denominator of the form 
factor, a = 2, is chosen to reproduce the Drell-Yan-West relation for large x. The mass difference 
M a — M s is motivated by the N — A mass difference (with group theoretical factors properly ac- 
counted for) and the values for A and M s reproduce the experimental value for the axial charge, 
9 A- , , 



parameter 


fixed by: 


a = 2 
M a -M s = 200 MeV 
M s = 600 MeV; A = 0.5 GeV 
N 


Drell-Yan-West relation: /"(x) ~ (1 — x) 3 
N — A mass difference 
9A = 1-26 

number sum rules 

Jdxmx) = 2;Jdxff(x) = l 



4 numerical results 

Having fixed the parameters of the model we can present our numerical results. 

• The transverse momen- 
tum dependent distribution 
fi(x,p T ) as obtained from 
the analytical expression for 
$ by tracing with 7 + and 
integration over dp~ . The 
dependence on p T is driven 
by the choice of the form 
factors Eq. (|26D and momen- 
tum conservation. 



N 2 5(p 2 -2P-p + M 2 - M 2 R ) 



2(2tt) 



\p 2 -K 



2 1 2a 




0.25 1 
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Integrating fi(x,p T ) over d 2 p T we obtain the usual distribution fi(x). The values for the 
first moments Jq dxxfi(x) are 0.690 and 0.256 for the u and d quark, respectively. We 
compare our result with the parametrization from Gliick, Reya and Vogt |10| (at the low 
scale fM 2 LO = 0.23 GeV 2 )F|. 



xfj(x) , 


„ xf/x) , 
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/ X x Uy(x,ii L0 2 ) [GRV] 


/ \ - -- d - quark 0.6 ■ 


/ \- xdJx,ii L0 2 )[GRV] 


1 - 


/ \ OA- 
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s'' \ \ 02 ' 






0.5 1 x 


0.5 1 x 



Note that the lowest moments, J j\(x)dx are exactly the same (normalization condition), a 
fact not immediately apparent from the diagram, since we plot the combination x * f\ (x) . 
The position of the maxima is in fair agreement. Our distribution is narrower due to 
the non-inclusion of gluons and anti-quarks (which - if included in our model - would 
have a broadening effect). Thus, we refrain from fine-tuning parameters to obtain a closer 
agreement with GRV — and are satisfied with agreement in the gross features. 



The comparison of our result for g±(x) with parametrizations taken from the literature [|TT 
reveals a similar agreement in the gross features 



xg,(x) , 




x g/x) / 


xAuJx,il L0 2 ) [GRSV] 


1.5 


/ \ u - quark 


OA - 


xAdJx,ii L0 2 ) [GRSV] 


1 ■ 


/ \ - - d - quark 
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-0.2 - 


as / v 



Having gained some confidence in the model predictions - as well as insight in limitations 
in accuracy - we can make predictions for the less well known or unknown functions. 



6 Note that the model does not provide a scale dependence; but it is expected to describe physics at a low 
: 'hadronic" scale. Thus we compare to distributions found in the literature at scales as low as available. 



xh/x) a 



1.5 
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0.5 



-0.5 




u - quark 
d - quark 



Yet experimentally completely 
undetermined is the trans- 
verse spin distribution hi(x). 
Within our model we obtain 
an estimate for this function 
which is numerically close to 
the result for the helicity dis- 
tribution gi(x), but not identi- 
cal. 



We now turn to the 'twist 3' distribution functions. We display e(x) and the combination 
02 0*0 = 9t(x) —gi(x). The combination h 2 (x) = 2{hi{x) —hi(x)) is predicted in the model 
to be just 2 x g2(x). 
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From hermiticity and invariance under parity transformation relations between 'leading 
twist' and 'sub leading twist' functions can be found like [O, RJ 



d 



9t{x) 



h L (x) 



S 1 ), 



9i{x) + fc9iT{ x ) 



d 



ax 



(31) 
(32) 



involving p T -moments of distribution func- 
tions 

g$(x) = J d 2 p T (j^J 9it(x,Pt) 

which, as well, are predicted by the model. 

The p T -moment hfj^\x) is given by the re- 
lation ^ (x) - 



~9vt ( x ) within this model. 

Note the non-vanishing values of (and 
correspondingly of ) at x — 0. 
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The positivity constraints |<7i(sc)| < fi{ x ) an d < fi( x ) are trivially fulfilled in the model 

as well as the Soffer |T3| inequality, 2|/i 1 (x)| < (fi(x) + gi(x)). 

sum rule 



On the other hand, the non-zero values of g\x(x = 0) and h^ 1J (x = 0) indicate a small violation 
of the Burkhardt-Cottingham |TJ[ sum rule, Jq 1 dx gi{x) = 0, and the Burkardt [[HJ sum rule, 
Jq dx hi (x) = 0. Those violations turn out to be proportional to quark mass effects. Also, 
the Efremov-Teryaev-Leader |16| sum rule, Y, a ev e a Jo dx x (g1(x) + 2g%(x)) = 0, is violated. 
Probably, those violations are artefacts of the model and might be used to constrain future 
refinements. 



5 conclusions 

In the framework of a diquark spectator model we have obtained analytical expressions for 
quark-quark correlation functions defined as hadronic matrix elements of non-local operators. 
Distribution functions are given as Dirac projections of the correlators. Numerical results for 
leading and, in particular, sub-leading 'twist' functions have been reported and discussed. 
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